Valence bond crystal in a pyrochlore antiferromagnet with orbital degeneracy 
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We discuss the ground state of a pyrochlore lattice of threefold-orbitally-degenerate S — 1/2 
magnetic ions. We derive an effective spin-orbital Hamiltonian and show that the orbital degrees 
of freedom can modulate the spin exchange, removing the infinite spin-degeneracy characteristic 
of pyrochlore structures. The resulting state is a collection of spin-singlet dimers, with a residual 
degeneracy due to their relative orientation. This latter is lifted by a magneto-elastic interaction, 
induced in the spin-singlet phase-space, that forces a tetragonal distortion. Such a theory provides 
an explanation for the helical spin-singlet pattern observed in the B-spinel MgTi2C>4. 

PACS numbers: 75.10.Jm, 75.30.Et 
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Geometrically frustrated antiferromagnets have 
reached an increasing interest in the past decade 0. 
The reason is that their ground states are highly 
degenerate and can evolve in a variety of ways: they 
can remain liquid down to the lowest temperatures due 
to quantum effects 0, or lift their degeneracy via the 
order-out-of-disorder mechanism Q, or through a phase 
transition that lowers the local symmetry of the lattice 
4, 5]. In this Letter we want to point out and discuss 
another scenario, that can appear when magnetic ions 
of frustrated lattices possess also an orbital degeneracy. 
The physical behavior of such systems may be drastically 
different from that of pure spin models, as the occurrence 
of an orbital ordering can modulate the spin exchange, 
thus lifting the geometrical degeneracy of the underlying 
lattice. In the following we focus on a system with 
threefold-orbitally-degenerate S — 1/2 magnetic ions in 
a corner-sharing tetrahedral (pyrochlore) lattice. This 
model is suitable to describe d 1 -type transition-metal 
compounds, like the B-spinel MgTi 2 04. Here the 
magnetically active Ti 3+ -ions form a pyrochlore lattice 
and are characterized by one single electron in the 
threefold degenerate £2g-manifold. MgTi204 undergoes 
a metal-to-insulator transition on cooling below 260 
K, with an associated cubic-to-tetragonal lowering of 
the symmetry [(J. At the transition the magnetic 
susceptibility continously decreases and saturates, in 
the insulating phase, to a value which is anomalously 
small for spin 1/2 local moments: for this reason the 
insulating phase has been interpreted as a spin-singlet. 
Subsequent synchrotron and neutron powder diffraction 
experiments have revealed that the low-temperature 
crystal structure is made of alternating short and long 
Ti-Ti bonds forming a helix about the tetragonal c-axis 
0- These findings have suggested a removal of the 
pyrochlore degeneracy by a one-dimensional (ID) helical 
dimerization of the spin pattern, with spin-singlets 



located on short bonds. This phase can be regarded as 
a valence bond crystal (VBC) since the long-range order 
of spin-singlets (dimers) extends throughout the whole 
pyrochlore lattice. 

Here we describe the microscopic theory behind the 
stabilization of this VBC ground state. We argue that, 
remarkably, such a novel phase can be realized on a py- 
rochlore lattice because of orbital degeneracy, without in- 
voking any exotic interactions. The existence of orbitally- 
driven VBC had been suggested for cubic lattice of d 9 
compounds 8]. For d 2 compounds with frustrated lat- 
tices, the orbital order is shown to induce a spin-singlet 
ground state, for triangular lattices [|J, or a spin ordered 
one, for pyrochlore lattices [ljj. Yet, the peculiar case of 
d 1 spinel compounds leads to new results: the onset of 
an orbitally driven VBC state on the pyrochlore lattice. 

Effective Hamiltonian: We first derive a superexchange 
spin-orbital Hamiltonian a la Kugel-Khomskii (KK) [llj 
for threefold orbitally-degenerate dVions on a pyrochlore 
lattice. We assume that the low-temperature insulating 
phase of MgTi204 is of Mott-Hubbard type. We work 
in the cubic crystal class and look for possible instabil- 
ities towards symmetry reductions. Our parameters are 
the nearest-neighbor (NN) hopping term t, the Coulomb 
on-site repulsions U± (within the same orbital) and U2 
(among different orbitals), and the Hund's exchange, Jh- 
For t2 g wavefunctions the relation U± = U% + 2Jh holds 
due to rotational symmetry in real space. The orbital 
occupancies of ti g orbitals, n Q( 3 (a, ft = x,y,z), are ex- 
pressed in terms of the pseudospin f = 1, with the cor- 
respondence: t z — —1 — > \yz), t z — — * \xy), and 
t z = 1 — > \xz). At first, we consider only the leading 
part of the hopping term, due to the largest dda element, 
and discuss later the effects of smaller contributions (e.g. 
ddit). The dda overlap in a(3 plane connects only the cor- 
responding orbitals of the same a(3 type. Thus, the total 
number of electrons in each orbital state is a conserved 
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quantity and the orbital part of the effective Hamiltonian 
H e g is Ising-like: 

fleff = -Jl^[^i- Sj + 3/4] Oij (1) 

(v) 

+J 2 ]T [3 ■ Sj - 1/4] O - + J 3 ^ • ^ - 1/4] Oij 

where the sum is restricted to the NN sites. Intro- 
ducing the projectors on the orbital states of site i, 
Pi,xz = §Ti Z (l + T iz ), P i<xy = (1 - T iz )(l + t iz ) and 
Pi,yz = — f r iz(l — r iz): the orbital contributions along 
the bond ij in a/3-plane is given by Oy = Pj ia/ g(l — 

Pj,ap) + -fj,a/3(l - Pi.a/s) and Oy = Pi, a pPj, a p. The 
first and second terms in H c s describe the ferromag- 
netic (FM) Ji = t 2 /(U~2 — Jh) and the antiferromagnetic 
(AFM) J 2 — t 2 /(U 2 + Jh) interactions, respectively, and 
are active only when the two sites involved are occu- 
pied by different orbitals. The last term is AFM, with 
J 3 = \t 2 [2/{U 2 + Jh) + l/(U 2 + 4J H )] , and is non-zero 
only when the two sites have the same orbital occupancy. 
At this point it is useful to have an idea of the energy 
scales that play a role in the Hamiltonian (JTJ. We es- 
timate t = t„ ~ 0.32 eV, J H ^ 0.64 eV and U 2 ^ 4.1 
eV H3. Thus 77 = J H /U 2 ~ 0.15 < 1 and, just in or- 
der to present the results in a more transparent form, 
we expand the exchange energies around r] = 0. We get 
Ji ~ J(l + rj), J 2 ~ J(l - rj) and J 3 ~ 4 J(l - 2rj) where 
J = t 2 /U2 — 25 meV represents the overall energy scale. 
In the following we measure all energies in units of J. 

The main aspect of H c r is that, due to cider-character 
of the hopping terms, only some orbital configurations 
contribute to the energy: every bond ij in the a/3 plane 
has zero energy gain unless at least one of the two sites i 
and j has an occupied a/3 orbital. The strength, as well 
as the sign, of spin-exchange energy associated with two 
NN sites i and j depends on their orbital occupations and 
the direction of the ij bond. The strongest bond in the 
generic a/3-plane is characterized by both sites with a/3 
occupancy: we shall call it bo- Its exchange interaction 
is AFM and its spin Hamiltonian is given by: 

H bo = -l + 2r 1 + 4[l-2r ) }S i -S j . (2) 

When the two sites of bond ij in a/3-plane are occupied 
by one 07 and one a/3 orbitals, 7 7^ /3 (bond bi), one gets 
a weak FM interaction: 

H bl = -1 - v /2 - 2r)Si ■ § s . (3) 

Finally, the two sites of bond ij in a/3-plane can be oc- 
cupied by one 07 and one /?7 orbitals (bond 62), or, by 
two 07 (or (3-y) orbitals (bond 63). These bonds are non- 
interacting, as far as only dda overlap is considered. 

Single tetrahedron: In one tetrahedron there are basi- 
cally three possible orbital configurations to be consid- 
ered (see Fig. [IJ: A) all four sites have the same orbital 
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FIG. 1: Orbital and bond arrangements on a tetrahedron for 
cases ^4), B) and C). 

occupancy (say xy) and thus only the two bonds in xy- 
plane (shown by solid lines in Fig. ^V) give a non-zero 
energy contribution; B) the two sites in one a/3-plane, e.g. 
xz 7 are both occupied by xz orbitals, while at least one of 
the two sites in the other xz-plane is occupied by xy or yz 
orbitals (Fig. ^3) ; C) no bonds ij in the plane a/3 of the 
tetrahedron is such as to have both sites occupied by the 
a/3-orbital (Fig. QJU). These three configurations are the 
bricks that allow to build the orbital pattern throughout 
the whole pyrochlore lattice. Because of the Ising-form of 
orbital interactions, in the following we can focus simply 
on these three cases, relying on the fact that configura- 
tions with a linear superposition of orbitals on each site 
must have a higher energy. We shall do only one excep- 
tion to study a case with a particular physical meaning, 
i.e., that of a "cubic" symmetry, where each site is occu- 
pied by a linear superposition with equal weights of the 
three orbitals -j^{\xy) + \xz) + \yz)], (case D). 

Pyrochlore lattice: Here we consider possible coverings 
of the pyrochlore lattice through the various tetrahedra. 

A) Heisenberg chains - When all tetrahedra of py- 
rochlore lattice are of type A (ferro-orbital ordering) then 
the effective Hamiltonian JJJ can be mapped into a set 
of one dimensional decoupled Heisenberg chains. If, for 
example, all occupied orbitals are of xy-type, all chains 
in (1, ±1, 0) cubic directions (see Fig. QJV) are decoupled. 
The only interactions are due to AFM 60-bonds described 
by the spin Hamiltonian J2J). Thus, the ground-state en- 
ergy per site can be evaluated exactly by using the results 
for an Heisenberg ID-chain 0: E A = -2.77(1 - 2rj). 

B) Dimer phase - This state is made of i?-type tetra- 
hedra. We distinguish three types of such tetrahedra. All 
three are characterized by one strong bo bond, in a a/3 
plane. The other two ions in the opposite a/3-plane can 
either be occupied by two (/3j) orbitals, forming a 63 
bond, (case B\, the one shown in Fig. HJ3), or by one 
07 and one a/3 orbitals, linked in a b\ bond, (case B 2 ), 
or, finally, by one orj and one /?7 orbitals, forming a b 2 
bond, (case -B3). Since b 2 and 63 bonds do no contribute 
to the energy, this latter depends only on the number, 
ni, and n^, of bo and b\ bonds in the unit cell. As all 
three Bi configurations are characterized by = 1 and 
ni >1 =2, all possible coverings of pyrochlore lattice by Bi 
tetrahedra have the same energy, even if (n(, 2 , ) are 
different for three Bi tetrahedra ((2,1) for B\\ (1,2) for 
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Bi\ (3,0) for B 3 ). When pyrochlore lattice is covered by 
-Bj-types tetrahedra (two possible coverings are shown in 
Fig. |2J), then each spin is engaged in one strong AFM 
b bond Eq. J2J), and two weak FM ^-bonds Eq. |J3J. 
Such coverings form a degenerate manifold and the corre- 
sponding energy can be calculated as follows. In the limit 
X] — *■ 0, the spin-only Hamiltonian can be solved exactly, 
as it can be decomposed into a sum of spin-uncoupled bo 
bonds. In this case the energy minimum is reached when 
the Heisenberg term of the &o - bond is the lowest, i.e., for 
a pure quantum spin-singlet (Si ■ Sj — —3/4). Remark- 
ably, such spin-singlet (dimer) states, in the limit 77 — > 0, 
are also exact eigenstates of the full Hamiltonian (JJJ . As 
i) « 1, the dimer state is stable against the weak FM 
interdimer interaction. In this case the magnetic contri- 
bution along the FM fei-bond is zero ((Si ■ Sj) = for i 
and j belonging to different dimers) and we are led to an 
energy per site given by: Eb = Et /2 + Et 1 = — 3 + |??. 
Here E(, 0(1) is the energy of the bond &o(i)- 

C) FM order: Consider the state where all tetrahedra 
are of type C (see Fig. \J]p). There are four interacting 
b\ bonds and two noninteracting bonds (62 and 63) per 
tetrahedron. All non zero spin-exchanges are FM given 
by Eq. (0). The ground state for this case is, thus, 
ferromagnetic with an energy Eq = 2Ef )1 = — 2(1 + 77). 

D) Frustrated AFM: The realization of this phase 
restores the full pyrochlore lattice symmetry, thus de- 
scribing an ideal cubic phase. All bonds are equiva- 
lent and by averaging Eq. over the orbital config- 
urations on neighboring sites i and j, we obtain the 
spin Heisenberg Hamiltonian on the pyrochlore lattice: 
Hd = E ( y>(-5/9 + [4/9 - 1677/9]^ • S,). The system 
is thus highly frustrated and its ground state is a spin- 
liquid |3, whose energy per site is: Ed ~ —1.89 + 0.89?7. 
Here we have used the ground-state energy estimated 
(l/N) Si ■ Sj - -0-5 on the pyrochlore lattice Q. 

E) Mixed A and C configuration: It is possible 
to cover the pyrochlore lattice also by means of a mixed 
configuration with A and C types tetrahedra. It can be 
visualized from Fig. |5Jd if, e.g., the four labeled tetrahe- 
dra would be of A-type, with two spin-singlets on strong 
bo bonds and three different orbital occupancy, and the 
four unlabeled tetrahedra of C-type with no spin-singlets. 
This dimer configuration is degenerate with the B-phase 
as far as only dda overlap is considered, as, in average, 
rib — 1 and n\ >1 = 2. Yet, in this case rib 2 — 1/2 and 
rit, 3 =5/2 and the degeneracy is removed by ddir overlap 
in favor of the B-phase (for which n^ 2 > 1), as the energy 
gain of the bi bond is —t^/[At 2 } (in units of J), while that 
of the b 3 bond is -t 2 n /[8t 2 }. 

Ground state manifold: On the basis of the previous 
energy considerations, a simple phase diagram can be 
derived, in terms of 77, the only free parameter available. 
For 37 = the lowest ground-state energy is that of phase 
B. With increasing 7/ we find only one phase transition at 
i] c = 2/11 ~ 0.18, between dimer-phase B and FM phase 




FIG. 2: Two different coverings of the unit cubic cell through 
dimers, phase B. The same notations as in Fig. 0are used. 
Locations of singlets are represented by thick links. Differ- 
ent numbers correspond to inequivalent tetrahedra. (a) The 
experimental phase of MgTi2 04: the helical dimerization pat- 
tern (indicated by arrows) is formed by alternating short 60 
and long 63 bonds. All tetrahedra are in phase Bi; (b) a differ- 
ent covering of the cubic cell through all Bi tetrahedra. The 
inclusion of the magneto-elastic coupling pushes these states 
higher in energy. 

C. As rj c is above our estimated value of 77 ~ 0.15, we can 
conclude that the ground state of MgTi204 is described 
by the phase B and is characterized by a frozen pattern 
of spin-singlets throughout the whole pyrochlore lattice 
that removes the original spin degeneracy. Nonetheless, 
there is still a remaining degeneracy to be lifted. It is re- 
lated to the freedom in the choice of the two orbitals on 
the tetrahedron bond opposite to the one of the singlet. 
Different choices of these orbitals give rise to inequivalent 
covering patterns of the pyrochlore lattice with one dimer 
per tetrahedron (see Fig. |2J). This degeneracy is given by 
the number of such dimer coverings and the correspond- 
ing number of states can be estimated to grow with the 

system size as Af ~ 3 Nt = VS* Here N T = N/2 is 
the number of tetrahedra and we have ignored the contri- 
butions coming from closed loops (hexagons) on the py- 
rochlore lattice. This ground state manifold is different 
from a resonating valence bond state, since each dimer 
covering is frozen in an exact eigenstate of the Hamil- 
tonian Q for 77 = 0. For finite 77 the different dimer 
patterns are not connected by Hamiltonian: the bond 
corresponding to the dimer in each tetrahedron is fixed, 
being determined by orbital pattern and orbital degrees 
of freedom are Ising-like variables. Thus a tunneling be- 
tween different dimer states cannot occur. 

Lifting of degeneracy: The above discussed degeneracy 
can not be removed within the KK-type model, not even 
introducing ddn and ddd overlaps. The reason is related 
to the fact that the energy gain depends only on the total 
number of each type of bond (n.& , n^, nb 2 , n& 3 ) in the 
unit cell and, in order to fill the whole crystal with a 
periodicity not lower than the one of the primitive cubic 
cell shown in Fig. [21 the average number of bonds ri6 i per 
tetrahedron is the same, whichever of the three building 
blocks Bi, B 2 , or B 3 is used. It is given by nb = 1, 
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rib 1 = 2, ri}y 2 = 2, ni, 3 — 1, and it corresponds to the value 
of case B\ , that is the only one that allows to cover the 
whole cubic cell without mixing to other configurations 
(see Fig. Associated to the i?i-phase, we have the 
minimal cell enlarging (doubling instead of quadrupling) , 
and the maximal space subgroup (P4i2i2) of the original 
face-centered Fd3m cubic cell. 

In order to substantiate these geometrical considera- 
tions, we need to find the physical mechanism that re- 
moves the B- manifold degeneracy in favor of the in- 
state. From the above discussion it follows that only 
correlations between bonds can lift it. These correlations 
naturally appear if the magneto-elastic contribution to 
the energy is considered, as the orbitally-driven modu- 
lations of the spin exchange energies distort the under- 
lying lattice through the spin-Peierls mechanism. In the 
degenerate B-manifold, every tetrahedron is character- 
ized by a strong exchange on the bond bo where the sin- 
glet is located. A reduction of the bond length enlarges 
the energy gain, because of the increase in the dda over- 
lap. This selects the triplet-T deformation mode from 
the irreducible representations of the tetrahedron group, 
which is the only one that singles out one shorter bond 
[lif . This mode generates a tetragonal distortion of the 
tetrahedron, with short and long bonds located opposite 
to each other and four intermediate bonds. Due to this 
mechanism, the position of the two &i-bonds in the tetra- 
hedron is uniquely determined: in order to maximize the 
superexchange energy gain by keeping the highest value 
for t a , the intermediate-strength &i-bonds are not allowed 
to lie on the long bond opposite to the singlet bo- The 
elongation of the weak bonds of 63 type is energetically 
more favorable. It is possible to check that the only pos- 
sibility to have such a constraint for the whole cell is 
realized for the state B\ (Fig. On the contrary, 

both cases B2 and B3 (Fig. Eb), do not allow to cover 
the cell without at least one 61-bond lying opposite to the 
singlet edge (e.g. tetrahedron 4 in Fig. 2b), thus with 
an extra-energy cost. Hence, the energy is minimized 
when all tetrahedra are of B\ kind, with a T-type tetrag- 
onal distortion. In this state all dimers are condensed in 
the ordered helical pattern shown in Fig. [2b- and form a 
VBC. This dimerization pattern exactly reproduces the 
one observed in the insulating phase of MgTi204 0- The 
present theory also predicts a peculiar orbital ordering in 
the dimerized phase: a ferro-type order along the helices 
with antiferro-type order between them (see Fig. E^). 
This orbital ordering can undergo an experimental test 
through Ti K edge natural circular dichroism, which is 
sensitive to the chirality of t^g orbital order along the 
helix, when x-rays are shone along the helical axis. 

Within the present scenario two transition tempera- 
tures are expected. The highest T Cl , determined by the 
exchange coupling within the singlet (T Cl ~ J3 ~ 4J ~ 
1000A"), corresponds to the transition from a paramag- 
net to a spin gap (dimer) state, with a ferro-orbital order 



on each dimer. This state can be regarded as a weakly- 
interacting gas of dimers and is highly degenerate with 
respect to the dimers orientation. The lowest transition 
temperature, given by the magneto-elastic coupling, lifts 
the degeneracy through the spin-Peierls distortion. At 
this temperature dimers condense and form the VBC 
shown in FigEJi. The entropy involved in this transi- 
tion is estimated to be ~ In [Af] /N = ln\/3. In the case 
of MgTi2C>4, it is known |y| that, with increasing temper- 
ature, this compound goes from an insulating to a metal- 
lic phase at T C2 ~ 260A. The transition to the metallic 
state rules out the possibility of a high-temperature spin- 
singlet state with disordered dimers and does not allow 
to evaluate the order of magnitude of the magneto-elastic 
coupling. We can just estimate its lower limit as about 
l/4(~ T C2 /T Cl ) of the singlet-binding energy. 

In summary, we have derived a mechanism that allows 
to lift the geometrical degeneracy of a pyrochlore lat- 
tice of threefold-orbitally-degenerate t 2g , spin- 1/2 mag- 
netic ions like Ti 3+ . We have singled out two relevant 
energy scales that govern its behavior: the main one is 
the superexchange energy gain, that drives the system 
into a spin-singlet dimer state with peculiar orbital pat- 
tern. The residual orientational degeneracy is then lifted 
through the magneto-elastic interaction that optimizes 
the previous energy gain by distorting the bonds in the 
suitable directions and leading to a tetragonal distortion. 
This generates a condensate of dimers in a VBC state, 
forming ID dimerized helical chains running around the 
tetragonal c-axis, the one actually observed in MgTi2C>4. 
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